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Abstract 

Application of the density functional formalism to the canonical en- 
semble is of practical interest in cases where there is a marked difference 
between, say, the canonical and the grand canonical ensemble (cavities or 
pores). An important role is played by the necessary modification of the 
famous Ornstein-Zernike relation between pair correlation and direct cor- 
relation function, as the former is no longer invertible in a strict sense in 
(finite) canonical ensembles. Here we approach the problem from a differ- 
ent direction which may complement the density functional approach. In 
particular, we develop rigorous canonical ensemble versions of the LMBW 
equations, relating density gradient and exterior potential in the presence 
of explicit (singular) containing potentials. This is accomplished with the 
help of integral operator and Hilbert space methods, yielding among other 
things representations of the direct correlation function on certain sub- 
spaces. The results are particularly noteworthy and transparent if the 
segregating potential is a linear (gravitational) one. In that case the mod- 
ifications in the LMBW equations can be expressed as pure, seemingly 
non-local integrations over the container boundaries. 



1 Introduction 



Recently there has been renewed interest in the the apphcation of the density 
functional method to the canonical ensemble (CE). One of the more practical 
reasons is the investigation of fluids in cavities or pores that are so small that 
the results depend on the choice of the ensemble (see, e.g., P or, for a short 
review, j2j). One of the methods that have turned out to be quite useful for the 
description of finite systems consists of the so-called density functional theory. 

However, in its customary form this method requires the use of a grand canon- 
ical ensemble (GCE), i. e., an open system. Consequently this restriction has led 
to the (re) consideration of the question in how far the methods of density func- 
tional theory can be transferred to the case of closed systems that have to be 
described in the framework of a canonical ensemble (CE) (c/. [3j and the more 



It is the constraint of a fixed number of particles which causes certain technical 
problems. The reason is discussed in the following. We note that these problems 
have already been observed in the past for several times but, as remarked by 
Ramshaw [3], have occasionally been lost sight of or glossed over. 

The density functional method is based on the observation that for a finite 
system the external potential, v{r), and the one-particle density, p(r), are stand- 
ing in a one-to-one relation in a gas or fluid (see pj). We note however that there 
may arise problems as to this property if one passes to the thermodynamic limit 
and one is in a regime where the system develops long-range correlations as, for 
example, in the situation of phase transitions or spontaneous symmetry breaking. 

In finite systems, however, this result leads via the functional variation of 
p(r) with respect to v{r') and vice versa to an important relation between the 
two-particle correlation function, H{r,r') and its inverse, C(r, r'), the direct cor- 
relation function. This relation is frequently called the Ornstein- Zernike relation. 

It was already observed by Mermin [H], that this uniqueness relation cannot 
hold in the CE (in a finite volume), as the density remains the same under a 
constant shift of the external potential. The direct consequence of this fact is 
that 



in the CE in a finite volume, V. Put differently, if we regard H, C as the ker- 
nels of the corresponding integral operators, H, C on some function spaces, the 
constant functions are zero eigenfunctions of H. This implies that it cannot be 
straightforwardly inverted and hence the direct correlation function, C, is (or 
rather seems to be) ill-defined (see appendix). This is, in a nutshell, the problem 
of the Ornstein- Zernike relation in the CE. 

The above phenomenon has been observed in the literature for several times, 
but, nevertheless, it seems that it is by no means widely known (as has, for ex- 
ample, been also remarked in pi). We think it is the underlying psychological 



recent papers [H], [E!, [ZD- 
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reason why people usually prefer to use the grand canonical ensemble without be- 
ing really aware of the above result (on the surface it is the fact that in functional 
derivations one has to keep the particle number constant). 

There have been several attempts in the past to remedy this problem (see the 
papers 0, jE], [Z| already cited above). As far as we can see, all of them seem 
to employ the functional variation method and arrive at more or less different 
recipes. 

We entered this field from a slightly different direction. Our main interest 
centered around the long-standing and controversial issue of the behavior of the, 
say, liquid-gas interface when treated from practically first principles, that is, 
within the regime of microscopic statistical mechanics. An important role in 
this context is played by the famous LMBW equations and their corresponding 
hierarchies [S], ^0] (see the next sections). The underlying common assumption 
in both fields is the above mentioned one-to-one dependence between exterior 
field and one-particle density. 

In ^T] we discussed in quite some detail what can happen if there exist long- 
range correlations in the system or in the interface between coexisting phases 
(as, for example, for vanishing gravitational field in the liquid-gas interface). We 
note that one cannot avoid these problems by making the calculations in finite 
systems and passing to the thermodynamic limit afterwards. In that case one 
has to control the uniformity of convergence of various expressions in this limit 
which is typically lost. For more discussion in this direction see also |12J or jT^] 
and further references therein. 

The present paper is organized as follows: In section 2 the first LMBW equa- 
tion (the one involving the density- density correlation H) is dealt with in the 
framework of the canonical ensemble. In particular, the exceptional significance 
of boundary terms brought about by the necessary confining part of the exterior 
potential is exhibited. Moreover, for the special case where the non-confining 
(segregating) part of the potential is chosen to be linear, it is shown that only 
boundary integrations survive on the r. h. s of the LMBW equation which cor- 
roborates the important role of long-ranged correlations and the delicacy of in- 
terchanging limits. 

In section 3 we derive a suitable generalization of the Ornstein-Zernike relation 
for the canonical ensemble case. In contrast to previous methods discussed in 
the literature, we refrain from explicit recourses to density functional theory. 
Extensive use is made instead of Hilbert space methods, the main point being 
the restriction of the non-invertible integral operator H to a subspace where it is 
again invertible. The correction terms for the Ornstein-Zernike relation appearing 
in our approach are then compared with those in the previous proposals 0,0) 

Our results allow in particular for an immediate derivation of a (generalized) 
version of the second LMBW equation (the one involving the direct correlation 
function C) for the canonical ensemble. This is done in section 4 and it turns 
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out that boundary terms have a great significance also for the case of the second 
LMBW equation. In fact, whenever the non-confining part of the exterior poten- 
tial is chosen to be linear, the integrations involved can be shown to completely 
reduce to boundary integrations again. 

The paper is then concluded with a summary. Furthermore, an appendix is 
devoted to the fixing of some elementary notions and notations. 

2 The First LMBW Equation in the Canonical 
Ensemble 

In the following we will study systems in a finite closed container, V, the walls be- 
ing explicitly incorporated via appropriate containing potentials. We now briefly 
dccribc how the containing potential is fed into the respective expressions. We 
split the total external potential into a segregating (bulk) potential and a contain- 
ing potential. 

v{r) ^ Vs{r) + Vc{r) (2) 

One of the technical advantages of this method is that one can perform the 
calculations as if one were in an inflnite system since Vc brings all quantities down 
to zero outside of V (see below). Furthermore, we can perform a rigorous scaling 
limit of Vc{r) near the walls, thus getting an exact expression for systems being 
contained by ideal hard walls. 

We start with a realistic wall potential, Vc, which is assumed to be only 
dependent on the coordinate normal to the wall and which is concentrated in 
a thin layer adjacent to the system boundary, dV, with Vc diverging smoothly 
to -|-oo at the boundary and being -|-oo outside V. For reasons of simplicity we 
illustrate the procedure with the help of a wall sitting at x — 0, the system being 
contained in the half-space x <0. 

Assumption 2.1 Vc{x) is assumed to be smooth for x < 0, Vc{x) being identically 

zero for x < —oq < and Vc{x) = +oo for x > 0. 

An inflnitely hard wall is represented by the scaling limit, A ^ oo, of 
v^{x) = Vc{X ■ x). It is instructive to control various quantities when one of 
the coordinates passes through the wall. Under these circumstances the expres- 
sions are dominated by the term exp — /3 • Vc{r), with Vc{r) — > oo at the wall. 
That is, we have for finite systems (where all the expressions can be rigorously 
controlled) , 

H{r, r') = f{r, r') exp {-[5vc{r')) (3) 

with f{r,r') continuous across the system boundary as long as r is sitting in the 
bulk. The same applies to p(r) for r passing through the wall. We have 

p(r) = /(r) exp (-/?Vc(r)) (4) 
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with /(r) again being continuous across the walL 

In contrast to most of the approaches to the topic described in the intro- 
duction, we want to choose the route via a rigorous derivation of the LMBW 
equations, which describe the relation between Vp(r) and Vf (r) by means of the 
correlation functions H{r,r') and C{r,r'). 

The LMBW hierarchies were described in jH], JU] with the help of density 
functional methods. It is not made entirely explicit what happens with boundary 
terms if one uses containing potentials and/or performs the thermodynamic limit, 
N,V ^ oo, in the presence of possible long-range correlations. Concerning this 
particular important point, we had an interesting discussion with R. Lovett some 
time ago ([E]; see also [IHI). He told us that [T3] was in part influenced by 
our earlier work in ^1] (see in particular section 5.3 of J3] to get a feeling of 
the epistemological problems underlying the need to perform a thermodynamic 
limit). 

In the following we show how to cope appropriately with the possible singular 
terms which may emerge in this limit. The important point is that, in order 
to enclose a finite system in a container (unless one uses for example periodic 
boundary conditions, which are a little bit artificial as long as one does not per- 
form a thermodynamic limit in the end), one should not gloss over the necessary 
bounding terms, which, as we showed rigorously, may make an effect even in this 
limit. 

In [TJj we started from the so-called Kubo-Martin-Schwinger (KMS) property 
of equilibrium systems which holds in finite as well as in infinite systems. We 
showed that both the BBGKY hierarchy and the first LMBW hierarchy emerge 
rigorously as the two opposite ends of one and the same scaling procedure, as long 
as certain cluster properties are fulfilled. If these cluster properties are violated, 
extra boundary terms may survive with possible severe consequences for concepts 
like the Triezenberg-Zwanzig formula for the surface tension. 

As a particular result of jllj it turns out that for finite systems held together 
by an explicit containing potential of the sort described above, the first LMBW 
equation rigorously holds in the canonical ensemble case (provided the containing 
potential is correctly dealt with). That is, we have for the first member of the 
hierarchy: 



with V = Vs + Vc- What interests us here is the term on the right side including 
the containing potential Vc- It reads 



We can simulate ideal hard walls by the method described above. Using 
again as illustration a wall, sitting at x = and Vc being supported in (— ao,0) 
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with Vc — oo ioT X > 0, we can rigorously perform the scahng process with 
v^{x) :— Vc{Xx) and get (for the wall described above): 

-13 J H{r,r')d:cV^{r')dr' = -13 J f{r,r') exp{-f3v^{x)) d^v^{x) dxclr'^"^-^^ (7) 

{d denoting the space dimension). 

Making a variable transformation, x' :— Xx, we get 

/(r, r'; x'/X)d^, exp{-Pvc{x'))dx'dr'^'^-^^ (8) 

where the meaning of the abbreviations should be clear. In the limit A — > oo we 
hence get for the x'-integration (with Vc{—ao) — , Vc{0) — +oo): 

J f{r, r'; x = 0) (^J exp {-(3vc{x))dx^ dr'^"^'^^ ^- J f{r, r'; x = 0) dr'^'^-^^ 

(9) 

Note that for a realistic (i.e., unsealed) containing potential, H{r,r') — > 
strongly for r' — > dV from the interior of y . In the case of an ideal hard wall 
(scaling limit) we have in general 

H{r,r') ^ H+{r,r) ^0 (10) 

for r' dV with the plus sign denoting the limit taken from the interior of V. 
For r' lying outside of V, H{r, r') is of course zero in both cases with usually 
having a jump-discontinuity. 

Conclusion 2.2 In case of an ideal hard wall the LMBW equation reads 

Vp{r) ^ -P f H{r,r')Vvs{r')dr' - [ H+{r,r')dd' (11) 

Jv JdV 

with do' the outwardly oriented (vector-valued) surface element. Note that 
H^{r,r') now depends only on the segregation potential Vg. 

Remark: We ignore the situation near possible points of the containing surface, 
dV, where the surface normal may change discontinuously. These points do not 
contribute in the scaling limit as they are of measure zero and can be avoided 
anyhow by making the surface smooth. 

Frequently the external gravitational potential is represented by a linear term, 
Vs{r) :— mgz, viz. 

Vvsij-) = m g Cz (12) 

[cz being the unit vector in z-direction). While for infinite systems one has to 
take into account the unboundedness from below of such a linear term, everything 
is well defined for a finite volume. 

Our above version of the first LMBW equation in the canonical ensemble has 
the following, as we think, interesting consequence. 
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Observation 2.3 As a result of the exact vanishing of Jy H{r,r') dr' and the 
assumed linearity of Vs, we have 

Vp(r) = -13 j H{r, r') Vv,{r') dr' (13) 

or 

Vp(r) = - [ H+{r,r')dd' (14) 

JdV 

Conclusion 2.4 As the containing potential is assumed to have its support re- 
stricted to a thin layer adjacent to the walls of the container, H or must 
include macroscopic (long-range) density correlations between the infinitesimal 
neighborhood of the container walls and the interior of the vessel as long as 
Vp(r) ^ 0. 

In this phenomenon was attributed to finite size effects in the canonical 
ensemble. A function like p^"^^ (r, r') with r' near or at the wall can be related 
to the conditional probability for finding a particle at position r in the bulk, say, 
under the assumption that another particle is sitting at a fixed position, r', near 
the wall. For free systems or short range correlations this quantity is expected 
to be (approximately) related to the two-particle distribution function for a (A^- 
l)-particle system within V. 

The consequence for the two-particle correlation function, H{r,r'), is that for 
large separation (r in the bulk, r' near the wall), z. e., |r — r'| — oo in a necessarily 
vague (or, rather, pragmatic) sense in a large but finite system, it does not go 
to exactly zero in a canonical ensemble but a finite size effect term of the order 
0{N~^) survives. 

These observations therefore furnish us with an exceedingly transparent ex- 
ample for which we can scrutinize some of the problems posed by Lovett and 
Buff in section 5.3 of In particular, it can now be clearly seen that there 
are indeed cases for which the structure at remote boundaries determines the 
local state of a system. The main point to be stated here is that this feature can 
nevertheless be brought into accord with the option to perform a thermodynamic 
limit. This is due to the extreme nonuniformity of the latter. For every sin- 
gle finite system with A^= const, the density gradient is expressed exclusively in 
terms of boundary integrations without explicit recourse to the segregating part 
of the potential due to the vanishing of J H{r,r') dr'. On the other hand, after 
having performed the thermodynamic limit, the information previously encoded 
in the boundary terms seems to have become unaccessible. However, it is by no 
means lost! It has merely been shifted to the segregating potential term which 
ceases to vanish in the N= oo case as the constant functions are no longer zero 
eigenfunctions of H there. 

We already emphasized in ^T] the important point that the interchange of 
limits in the derivation of the LMBW equation may be a delicate matter, as 
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the various limits might behave in a highly nonuniform way. Here we have seen 
that in the canonical ensemble such nonuniform behavior is even indispensable 
in order to have meaningful thermodynamic hmit states! 

Remark: As was already mentioned, a linear segregating potential is un- 
bounded from below, and it is therefore necessary to retain parts of the boundary 
walls in order to keep the system stable as N, oo. This does however not in- 
validate the above considerations. 

For the ideal gas these effects can be easily calculated. Without an external 
potential Vs{r) we have 

p^^\r, r') - p{r)p{r') = V'^iN^ - N) - V'^N^ = -{N/V) ■ ~ N'^ (15) 

With Vs{r) non-zero (and r 7^ r') we get: 

H{r,r') = p^'^\r,r')-p{r)p{r') = -N eK^{-l3{v{r)+v{r')))-{ J eicp{-Pv{r))dr)-'^ 

(16) 

The details of the scaling with the size of the system can be inferred from the 
following set up. We take for example a rectangular parallelepiped with area of 
the base equal to A, height equal to h and with the base located ai z = 0, the 
exterior potential being only 2;-dependent. We arrive at: 

H{r, r') = -N ■ A-'^{ J ex.p{-Pv{z))dz)-^ ■ exp(-/3(^;(r) + v{r'))) (17) 

and in the particular case v{z) — mgz this is equal to 

-A^ • A-'^ ■ {13 ■ mgf ■ (1 - e'^^'^^^'^)-'^ ■ exp(-/?(m^2 + mgz')) (18) 

Inserting the corresponding terms in the expression for i/(r, r'), r ^ r', in the 
grand canonical ensemble we get on the other hand for the ideal gas: 

p(2)(r,rO -p(r)p(r') = (19) 

If we take the limit /i — > 00 with N/A kept fixed, we end up with the baro- 
metric formula: 

H{z, 0) = (N/A) ■ A-^ ■ {(3mg)^ exp{-/3mgz) (20) 
That is, H scales as A"^. The LMBW equation leads to 

d/dz{p{z)) = (N/A) ■ {(3mgf ■ exp{-i3mgz) (21) 

or 

p{z) = {N/A) ■ {pmg) ■ exp{-pmgz) (22) 



8 



Conclusion 2.5 The local law between density and exterior potential is in the 
LMBW equation hidden in the non-local integration over the container walls. It 
can however be shown that in the case of the ideal gas this reduces basically to an 
application of Gauss 's divergence theorem. 

Remark: Even in the case of a zero external field and ideal hard walls the in- 
dividual contributions of the walls on Vp(r) in the bulk are non-vanishing and 
of mean-field type. In a finite volume and the ideal gas, they however exactly 
compensate each other such that we have in fact Vp(r) = as is the case by 
calculating p(r) directly. Once again we see from this that care has to be taken 
when performing the thermodynamic limit as the contributions do not decay and 
are independent of the distance of the respective walls! 

3 The Direct Correlation Function in the Can- 
onical Ensemble 

The ordinary Ornstein-Zernike relation 

J C{r, r") H{r", r') dr" = 5{r - r') (23) 

can be rewritten in integral operator notation as 

C H = 1 (24) 
where the integral operators C and H are defined as follows: 

(C/)(r) = J C{r,r')f{r')dr'; (H/)(r) = j H{r,r') f{r') dr' . (25) 

That is, C is nothing but the (left) operator inverse of H in some properly chosen 
linear space. 

Some remarks are here perhaps in order concerning the reason to introduce 
Hilbert spaces in the following. In contrast to, say, quantum theory, there is no 
inner theoretical reason to employ the Hilbert space formalism. The operator 
H lives on a certain space of functions. As the corresponding integral kernel, 
H{r,r') strongly decreases like exp{— l3vc{r')) when r' approaches the container 
walls, even singular functions like Vf (r') are in the domain of definition of H (see 
below). On the other hand, typically C(r, r') diverges with the same degree when 
r' approaches the container walls. Hence it is only defined on a smaller space of 
functions than for example L'^{V). We conclude that H naturally lives on a 
function space larger than L'^{V), while C lives on a smaller subspace of L'^{V). 
In any case, C is defined on the range of H which is dense in L'^{V) whenever 
H does not possess (proper) eigenfunctions with eigenvalue zero. However, the 
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validity of this latter property can only be expected if one works within the grand 
canonical ensemble. 

If one works in the canonical ensemble instead and again chooses = L'^{V) 
as Hilbert space, it is immediately obvious that H cannot be inverted there, 
since — as was already pointed out in the introduction (see also the appendix) — 
direct consequence of 

j H{r, r') dr' = (26) 

it possesses the constant functions as proper eigenfunctions with eigenvalue zero. 

However, due to the fact that we introduced a Hilbert space of functions as 
the domain of definition for H we now have the important technical notions of 
orthogonality and symmetric operators at our disposal, and we hence automati- 
cally know that eigenspaces are orthogonal and left invariant by the corresponding 
symmetric operator. In particular, these properties straightforwardly lead to the 
finding that a restricted form of invertibility for H is still present. To demon- 
strate this, let us carry through the just mentioned standard decomposition of the 
Hilbert space ^ into a direct sum of the eigenspace belonging to the eigenvalue 
zero and its orthogonal complement S^^: 

i3 = -^0 © i3o • (27) 

As H{r,r') = H{r',r), the integral operator H is symmetric in S). Therefore it 
leaves invariant the subspace 9)q and is invertible there. 

Under the premise that p[v] and p[v] differ whenever the external potentials 
V and V differ by more than a constant |H], it seems natural to assume that only 
the constant functions are annihilated by H, i. e., that is one-dimensional. In 
this case Sjq is given by 

y"/(r)dr = o|. (28) 

(The Hilbert spaces considered in this section consist of course of equivalence 
classes of functions rather than functions itself. For simplicity reasons we suppress 
this fact here.) 

A problem which we already addressed is that the Hilbert space = L'^{V) 
does not exhaust the full set of physically relevant functions on which one wants to 
define the operator H [cf. for example the LMBW equations). That is, in order 
to incorporate also these singular functions into the Hilbert space formalism, 
which is desirable on physical grounds, we have to modify the above ordinary 
Lebesgue measure. 

The mentioned divergencies can, for example, be compensated by the vanish- 
ing of the density p(r) which roughly behaves as exp{—j3 Vc{r)) in the vicinity of 
the container boundary. That is, one possibility consists of replacing the Hilbert 
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space L'^{V) with the Hilbert space L'^{V, p{r)dr) in which the measure is now 
chosen as p(r) dr instead of dr. While not being in L'^{V), the containing poten- 
tials and their derivatives are now naturally contained in L'^{V, p{r)dr). 

One should however emphasize that this choice is evidently not unique. Other 
options like, e. g., the replacement of p{r) with the corresponding ideal gas density 

= f 7'''7''';'' (29) 

J exp(— p v[r)) dr 

or with the quantity N{dp{r) / dN)^ (provided the latter expression is well-defined 
and positive) are possible. As we shall see below, this latter choice would in 
particular lead to results which are in some sense comparable with the treatment 
in 0, 0. 

In order to keep matters as general as possible, we therefore introduce the 
Hilbert space 

9) = L\V, p{r)dr) (30) 
where the measure p(r) dr fulfils the normalization condition 

p{r)dr = N. (31) 



Moreover, we require p(r) to be capable of compensating the singular behavior 
of Vf at the container boundaries. Thus, each of the above-mentioned quantities 
p, Pid, N{dp/dN)^ turns out to be a possible candidate for p. 

As a consequence, while H{r,r') is unique by its definition as a canonical 
expectation value of a certain microscopic observable and does not depend on 
the choice of the Hilbert space measure, the situation turns out to be quite 
different for C{r,r'). We think this observation is one of the virtues of our more 
functional-analytically oriented approach. 

An additional interesting feature of the Hilbert space Sj shows up if one deals 
with systems that still consist of a finite number of particles N but for which 
the volume is infinite. Such systems still make sense, e. g., in those cases where 
the confinement of particles comes about by the use of polynomial wells, i. e., 
potentials v{r) which diverge ~ \r\^ as |r| — > oo. 

The constant functions are no longer elements of in such a = oo case, but 
the integral operator H can still be applied to them and they are still annihilated. 
Therefore, the eigenvalue zero would be an improper eigenvalue lying outside the 
discrete part of the spectrum. 

In contrast to this, the constant functions are still contained in S^: 



p(r) dr = N<oo. (32) 
Therefore, zero is still an ordinary (proper) eigenvalue of H here. 
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The price one has to pay for this adaptation of the Hilbert space consists of 
the fact that the integral operator H is no longer symmetric, since 

</,H^>= J f{r){llg){r)p{r)dr = J f{r) H{r,r') g{r') p{r) dr' dr (33) 

is different from 

<H/,^> = j{nf){r')g{r')p{r')dr' = j f{r) H{r,r') g{r') p{r') dr' dr . (34) 

Among other things, this leads to the unpleasant feature that the subspace 
in the self-explanatory decomposition 

^ = ^00^^ (35) 

is not left invariant by H which means that a straightforward inversion is no 
longer possible there. 

However, one immediately notices that the modified integral kernel 

H'(ry)^^l^ (36) 

leads to a symmetric integral operator H', since 

< /, Ji'g>= J f{r) (H'g()(r) p{r) dr = J f{r)H{r,r')g{r') dr' dr 

= J f{r) H{r',r) g{r') dr' dr = J (H7)(r') g{r') p{r') dr' = < U'f,g > . (37) 

Moreover, the eigenfunctions of H and H' with eigenvalue zero coincide. Thus, 
as H' is symmetric, it leaves SJq invariant. Restriction to this subspace therefore 
makes an inversion of H' possible. 

As the elements of 9)o are annihilated by H', one therefore gets the following 
modification of the Ornstein-Zernike relation for the canonical ensemble case: 

C'H' = P^ (38) 

with P-*- denoting the projection operator onto the subspace JQq . 
Since 

(C'H7)(r) = J C'{r,r')H'{r',r")f{r")dr'dr" 

12 



it is immediately clear that the direct correlation function C{r,r') and its corre- 
sponding integral operator C have to be introduced via 

Thus, the modified Ornstein-Zernike relation for the canonical ensemble can also 
be rewritten as 

CH = P^. (41) 

Remark 3.1 While an inverse of H' does not exist outside a symmetric 
extension of the operator C can be defined on at least a dense subset of the full 
Hilbert space S), namely the linear combinations of elements from the range o/H' 
(which is at least dense in SJq ) and ioo- In particular, S^q is always left invariant 
by such an — of course non-unique — extension of C . 

If one again assumes that Jqo is one-dimensional, i. e., that H and H' only anni- 
hilate the constant functions, this implies 



C"(r,r')(ir' = const. (42) 
and hence we have the useful identity 

C(r, r')p(r') dr' = const. (43) 

Furthermore, i^o is then given by 

^o={/ <l,/>=y/(r)p(r)cir = o| (44) 

and P-*- is defined as follows: 

(PV)(r)=/(r)-^^^ = /(r)-l J f (r) ~p{r) dr . (45) 

The modification of the Ornstein-Zernike relation for the canonical ensemble can 
therefore be expressed as 

j C(r, r") H{r", r"')f{r"') dr" dr'" = /(r) j f{r) p{r) dr . (46) 

If one replaces /(r) with the delta function 6{r — r') concentrated at an arbitrary 
but fixed point r', the following final form emerges: 
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Conclusion 3.2 A modified Ornstein- Zernike relation in the canonical ensemble 
can be written down as follows: 



C{r, r") H{r", r') dr" = 6{r - r') - ^ . (47) 

Roughly speaking, our version of a modified Ornstein- Zernike relation differs 
from the one given in [EI by the appearance of the term p{r')/N instead of 
{dp{r')/dN),. 

To be more specific, let us compare our relations (jlSjl . (jUj), i. e., 

C(r, r") H{r", r') dr" + ^ = 6{r - r') , (48) 

C(r, r') dr' = const. (49) 

with the generalized Ornstein- Zernike relations that have been given as equations 
(52) and (53) in 0: 

C(r, r") H{r\ r') dr" + f = 5{r - r') , (50) 



ON 



V 



The difference of these latter relations to ours above is twofold. On the one hand 
it consists in the already mentioned replacement of p{r')/N with {dp{r')/dN)y. 
This difference reflects the freedom of choice for the measure of the Hilbert space 
Sj in our approach and would become inexistent with the particular definition 
p{r) dr = N{dp{r)/dN)y dr (provided the expression on the r. h. s. is a well- 
defined positive measure at all). On the other hand it turns out that the constant 
on the r. h. s. of ()5H1 is no longer arbitrary but has a definite value (correspond- 
ing to a definite extension of the operator C). Both differences can therefore be 
characterized as an implication of some freedom of choice in our approach which 
is not present in that of refs. |5j and 0. 

A brief remark is here in order concerning the results in (Hj, [7] and their 
relation to that of 0- We think that in both approaches the guiding idea 
is to mollify the canonical ensemble in order to make it more akin to the grand 
canonical case. This can be achieved in, on the surface, slightly different but 
nevertheless closely related ways. The strategy pursued in (5j, |j6j is to embed 
the canonical ensemble with fixed particle number into an enlarged parame- 
ter space, thus making a variation with respect to possible (see, however, the 
related short discussion in our summary at the end of the paper). Another pos- 
sibility is to make the system resemble an open system by truncating the weakly 
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decaying tails (coming from the constancy of A^) in the correlation functions as 
in jS] or jTj. Both strategies lead to closely interconnected results (c/. section 
3.4 of "G]) and correspond, as we have seen, to a particular choice in our own 
approach. 

However, nonwithstanding the fact that equations dSDI), ()51|1 are the natural 
outcome of the inversion formalism of the extended variable approach in |5, and 
it should be stressed that there is no obligation to make the correspond- 
ing choices for the measure and the constant also in our approach. Quite to the 
contrary, there is a certain amount of freedom which allows us to adopt some con- 
siderably simpler expressions. Furthermore, in the following we shall frequently 
deal with expressions where the operator C is solely applied to elements from 
the range of H'. In these cases we need not even worry about definite extensions 
of C. It is therefore obvious that — contrary to H{r,r') — the direct correlation 
function C (r, r') is by no means a uniquely defined quantity in our approach. 
This has the advantage that it is to some extent possible to simplify matters by 
using "tailor-made" definitions for C. 

To further clarify the contents of the present chapter, it might be useful to 
deal with the ideal g clS clS db relatively simple but thereby also fully tractable 
example. The expressions p(r) dr, Pid{r) dr, and N{dp{r) /dN)^ dr all coincide in 
this case and represent an obvious choice for our measure pdr. In the canonical 
ensemble we have 

Hid{r, r') = pid(r) S{r - r') - ^ pid(r) Pid(r') . (52) 

With p = we are therefore led to 

, n Hid{r,r') pid(r') 
Hid{r,r) = ^^-^ = 6{r-r) ^ (53) 

which means that i^o is indeed one- dimensional and that H-^ is simply given as 

Hid = P^. (54) 

After restriction to the subspace JQq , we thus get HJ^ = 1 and its inverse is 
trivially given as C(d = 1 there. This immediately leads to 

Cl^{r,r') = 6{r-r'). (55) 

For the time being, this expression is only valid on Sjq. However, its validity can 
immediately be extended to the full Hilbert space Jo (see remark EH}. Thus, in 
this case the direct correlation function reads 

C„(,.',^^^^^^^. (56) 

Pidin Pid{r) 
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Hence we ended up with exactly the same expression as for the grand canonical 
ensemble. Of course one has again to keep in mind here that the above extension 
of from 9)q to is by no means a unique one. So it should not come as a 
surprise that, e. g., the direct correlation function Cid{r,r') arrived at in equation 
(73) of P differs from our expression by a constant. 

As a further result we notice that the range of H(j consists of the full subspace 
SJq. This in turn allows for a definition of on the full Hilbert space Sj. Thus, 
the reservation "on at least a dense subset" made in remark 13.11 is unnecessary 
here. 

It is therefore interesting to deal with the question as to how far this desirable 
property is also transferable to the nonideal case. It will turn out in the following 
that a positive answer can be given for a rather large class of systems if we make 
the special choice p{r) = p{r). For a thorough discussion of the mathematics 
involved the reader may consult, e. g., ref. 

For p{r) = p{r) we have (see appendix): 

H'ir, r') = EIIjII = S{r - r') + h(r, r') pir') . (57) 
p(r) 

This corresponds to the natural decomposition 

H' = 1 + h' (58) 

with 

(h7) (r) = J h{r, r') f{r') p(r') dr' . (59) 
We now require h' to be a Hilbert-Schmidt operator on i^, i. e., 

j \h{r, r')\^ p{r) p{r') dr dr' < oo . (60) 

Under this proviso the Hilbert-Schmidt theorem tells us that h' is a completely 
continuous symmetric operator possessing a purely discrete spectrum with zero 
as its only possible limit (accumulation) point. Thus, the symmetric operator 
H' = 1+h' also has a purely discrete spectrum, but now with 1 as its only possible 
limit point. As an immediate consequence, the inverse of H' is defined on the full 
subspace iOg because the critical eigenvalue of H' is the eigenvalue zero which, 
according to our preceding discussion, can at most be finitely degenerated. The 
above result is also a particular example for the so-called "Fredholm alternative" 
which states that whenever k is a completely continuous operator on a Hilbert 
space, then either (l-fk)"^ is defined everywhere or l-|-k possesses an eigenvector 
with eigenvalue zero. 

It should finally be remarked that the simple form of the above natural de- 
composition (j58p of H' highly depends on the particular choice of g{r) dr as our 
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Hilbert space measure. For any other choice of the measure a nontrivial mul- 
tiphcation operator would show up instead of the 1, thus making an analogous 
reasoning as above more complicated. In particular, equation (j57j) had to be 
replaced with 



ZT// /X H{r, r') p{r) . 



(61) 

and it turns out that the further discussion would crucially depend on the prop- 
erties of the multiplication operator \/ pjp in the general case. 

Thus, again we see that matters can be kept simple due to the freedom of 
making the most suitable choice for the Hilbert space measure. 

4 The Second LMBW Equation in the Canoni- 
cal Ensemble 

Now we are ready to examine what expression replaces the second LMBW equa- 
tion (the one including the direct correlation function) when the canonical en- 
semble is used instead of the grand canonical one. 

As was already pointed out in section 2, the first LMBW equation 

(H Vt;)(r) = y E(r, r') Vv{r')dr' = Vp(r) (62) 

is still valid when the canonical ensemble is used. Due to the containing part of 
the potential, the derivatives of v{r) are strongly divergent in the vicinity of the 
container boundary. The replacement of Sj with S) therefore finds its justification 
here. Again putting Jq = L'^{V, p{r)dr), we get 

(H'Vt;)(r) = -l^. (63) 

The operator H' transforms elements of into elements of Sjq, and this is 
reflected in the fact that 

< 1,— >= /" Vp(r)rfr = (64) 
P J 

due to the vanishing of p(r) at the container boundary. 

The operator C can safely be applied onto both sides of the above LMBW 
equation and this leads to 

(C'H' Vt;)(r) = -i (c'^)(r) = (C Vp)(r) . (65) 
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Employing the modified Ornstein-Zernike relation C'H' = P-*-, we are thus im- 
mediately led to 

(C Vp) (r) = -p (P^Vv) (r) (66) 
or, if J5o is again assumed to be one-dimensional, to the following 

Conclusion 4.1 The second LMBW equation in the canonical ensemble can he 
written as 



j C{r, r') Vpir') dr' = -{3 (vv{r) - j P(r) Vf (r) dr I . 



(67) 



Therefore — in contrast to the unaltered first LMBW equation — this second 
one may differ from its grand canonical counterpart. The difference, however, 
ceases to exist for those choices of the measure where P-'-Vf = Vw, i. e., where 
the components of Vf are already elements of S^q. 

An example for this latter property is provided by the choice p(r) = Pidir) as 
equation immediately implies 

1 f 

Pid{r)Vv{r)dr = —— / Vpid{r)dr = (68) 



due to the vanishing of pid at the container boundary. 
There is, however, an even more remarkable example: 

Observation 4.2 Within the canonical ensemble we have 

j p{r)Vv{r)dr = 0. (69) 

Hence, for the choice p{r) = p{r) the components of'S/v{r) already lie in the or- 
thogonal subspace S^q. Since Vf contains both the segregating and the containing 
potential, we also have 

p{r) Vvs{r) dr = — j p{r) Vvc{r) dr. (70) 

Proof: The explicit inclusion of the containing potential greatly facilitates 
the calculations as we can perform all manipulations in the full M" (with the 
Boltzmann weight being zero for all configurations with the position of at least one 
particle lying outside the volume V). We assume that the interaction potential 
f/(ri, . . . ,rjv) between the N particles is translation invariant. By Va we denote 
the exterior potential translated by an arbitrary vector a G M": 

Va{r) = v{r + a) . (71) 
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Under these assumptions a simple variable transform shows that the canonical 
partition function QAr[t'] remains invariant under translations of v: 



QN[v]=QN[Va]. (72) 

The same holds for the logarithm and hence for the free energy. Differentiation 
with respect to a at a = immediately yields the desired result. q. e. d. 

The above result can also be shown in a slightly different way by employing the 
following equation which — again for the case of a translation invariant interaction 
potential — states that the time derivative of the momentum density must have 
zero expectation value in thermal equilibrium (c/., for example, §2 of [18j): 

Div7r(r) +p(r) Vt;(r) = 0. (73) 

Due to the appearance of the divergence of the pressure tensor 7r(r) in this 
equation, J p(r) Vf (r) dr can be reexpressed as a surface integral. Shifting the 
surface of integration beyond the container boundaries one again sees that this 
integral is in fact zero. 

The physical ramifications of this result are interesting. They show that 
the condition of equilibrium entails an a priori compensation of the global force 
exerted by the segregating potential on the system and the forces of constraint, 
provided by the container walls and acting on the nearby elements of the fluid 
or gas. This can most easily be seen in the case of the ideal gas (see our above 
calculations). 

Moreover, it is immediately seen that in those cases where well-behaved dif- 
ferentiation of / p(r) Vf (r) dr with respect to is possible, we also have 



dp{r] 
dN 



Vv{r)dr = 0. (74) 



Thus, whenever p(r) = N{dp{r) / dN).u is a possible choice, we have a further case 
for which P-'-Vf = Vf . 

In concluding this section we just give one little example for the implications 
which are brought about by the above results. Let the external potential as usual 
be divided into a segregating and a containing part. In particular, we choose the 
segregating part to be a linear (gravitational) potential : 

v{r) = Vs{r) + Vc{r) = mg-r + Vc{r) . (75) 

(Here and in the following g denotes a constant vector which can, e.g., be used 
to define the 2;-direction.) 
Then, due to 

Vfs(r) ~ / P(^) Vfs(r) dr = mg — — / p(r) dr = mg — mg = , (76) 
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the segregating part of v{r) drops out from the right hand side of the second 
LMBW equation and we are left with 



(C Vp)(r) = -P (vvc{r) j p{r) Vvc{r) dr^ . 



(77) 



(Strictly speaking, the segregating potential remains of course still present to 
some extent in those cases where we choose a p that implicitly or explicitly 
depends on it.) 

The containing potential Vc{r) can be chosen such that it identically vanishes 
for r outside an arbitrarily small vicinity of the container boundary dV. If r is 
sufficiently distant from this boundary, it therefore lies outside the support of Vc 
and the last equation reduces to 



(CVp)(r) = ^J p{r)Vv,{r)dr. (7J 



As in the grand canonical ensemble, where (CVp)(r) = —jSmg for r suf- 
ficiently distant from the boundary, (C Vp)(r) is still constant in the canonical 
case, but the value of this constant can be quite different in general. With special 
choices such as, e. g., 

it can even be made to vanish! 

Moreover, as we have already seen, 

^ J Po{r) Vv{r) dr^mg + ^ J po{r) Vvc{r) dr ^ (80) 

whenever po denotes any of the expressions p, pi^, N{dp/dN)v. This allows us to 
express the difference with the grand canonical case somewhat more explicitly: 

(C Vp)(r) = -Pmg + (p(r) - po(r)) Vv,{r) dr . (81) 

This leads to the interesting fact that the whole information about the value 
of (CVp)(r) in the canonical ensemble as well as about its possible difference 
with the corresponding grand canonical ensemble value —(3mg is contained in 
the values of certain functions on the support of Vvc, i- e., in an extremely small 
vicinity of the container boundary. 

This last statement can even be made somewhat more explicit by the obser- 
vation that one of course can also make the transition from smooth containing 
potentials to ideal hard walls here. Proceeding in an analogous way as in sec- 
tion 2, we are left with 

(CVp)(r) = l/ p^{r)do=-(3mg + ^ I {p+{r) - po + {r)) do (82) 

JdV JdV 

where the index once more denotes the limit taken from the interior of V . 
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5 Summary 



In the preceding sections we developed a method to estabhsh a kind of Ornstein- 
Zernike relation in the canonical ensemble which is complementary to the more 
common density functional approach. In the canonical ensemble the main prob- 
lem is that the two-particle correlation function, when viewed as an integral op- 
erator H, inevitably has a zero eigenvalue and thus cannot be straightforwardly 
inverted. Its inverse C, representing the direct correlation function, therefore 
seems to be ill-defined. 

Several suggestions have been made in the past to deal with or to circumvent 
this problem. In contrast to most of the treatments, wc emphasize an explicit 
operator theoretic point of view. We show that H can be inverted on some ap- 
propriate subspace of a certain function space being endowed with an adapted 
scalar product and being orthogonal to the constant functions which represent 
the eigenfunctions with eigenvalue zero. We introduce several intermediary func- 
tion spaces by means of which it becomes particularly transparent that (and 
moreover in what precise sense) C is non-unique on the full function space, the 
non-uniqueness stemming both from the chosen particular intermediary function 
space and the genuine non- uniqueness of the operator extension from a subspace 
to the full function space. 

We note that these intermediary spaces are practically unavoidable if one 
wants to apply H and C to their natural space of functions which comprise 
also singular functions like, e. g., Vfc(r), Vc denoting the containing potential 
confining the system to a finite container. These latter functions occur in the 
LMBW equations which represent a natural context for such relations and which 
we therefore treated as a straightforward application. To put it in a nutshell, 
while H is genuinely unique due to its explicit expression in statistical mechanics, 
this is not so for C. Different intermediary spaces and different extensions lead 
to slightly different expressions for the direct correlation function which however 
are all fulfilling a restricted Ornstein-Zernike relation for the canonical ensemble. 

A connection of our approach with the ones developed earlier in the realm 
of density functional theory can be made in those cases where the quantity 
{dp{r)/dN)v is well-defined and positive such that N{dp{r)/dN)y dr is a pos- 
sible choice for the measure of our Hilbert space. 

In this context it should perhaps be pointed out that the uniqueness or even 
well-definedness of quantities involving partial derivatives with respect to N is by 
no means self-evident. While the formal differentiation with respect to can be 
given a qualitative meaning in the realm of phenomenological thermodynamics 
{i. e., in some kind of thermodynamic limit where one can make use of the equiva- 
lence of ensembles) , we think the situation is much less clear whenever one wants 
to generate such expressions without transcending the realm of the canonical en- 
semble. This is even more so if one aims at the description of systems where the 
particle number is so small that the difference of ensembles becomes perceptible. 
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The main problem occurring here is that the particle number shows up in the 
dimension of the phase space volume and thus in the phase space measure as 
basically a discrete number. It is therefore not entirely evident how one should 
technically define unique continuous derivatives with respect to N without dis- 
missing the confinement to the canonical ensemble case. This problem was also 
observed by Ray JH] who even claims that there is no useful way to carry out 
such differentiations in general. 

One of the virtues of our approach therefore consists in the possibility to 
circumvent the above problem by choosing simpler measures such as p(r) dr or 
Pid(r) dr for our Hilbert space. Furthermore, differentiation with respect to N 
also plays a decisive role for the constant on the r. h. s. of equation (f^T|) . However, 
the non-uniqueness of the direct correlation function in our approach also enables 
us to alter this constant. It is therefore again possible to avoid differentiation 
with respect to N here. 

As a further advantageous feature of the mentioned non-uniquenesses we fi- 
nally note that the corresponding freedom of choice quite often allows for certain 
simplifications of procedures. 



Appendix 

For the convenience of the reader and in order to fix the notation we provide a 
short list of general relations. The pair correlation functions, H{r,r'), h{r,r'), 
are defined as 

H{r, r') := p^^) (^r, r') + p(r) 6{r - r') - p{r) p(r') , (83) 



h{r,r):= r^TT\ =-9{r,r)-l 84 

p{r)p{r') 

with the ordinary molecular distribution functions p(r), p^^^'^r, r') (see for example 
|2()j). The physical meaning of H{r,r') can be inferred from the relation 

H{r,r') = {6n{r)6n{r')) (85) 

with (...) denoting the ensemble expectation value and 6n{r) being the true 
microscopic observable 

6n{r) := n(r) — {n{r)) = n{r) — p{r) , n(r) := 6{r — rj) . (86) 

i 

This identity exhibits the fluctuation content of H{r,r'). 

The Ornstein- Zernike relation which in its general form can only hold in the 
grand canonical ensemble (see below) defines the direct correlation function and 
reads 

^ C{r, r') H{r', r") dr' = S{r - r") (87) 
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or with 

C(r, r') := 5{r - r')/p{r) - c(r, r') (88) 

we have 

Hr, /') - ^ / M, .>(0c(.', ^ / r'Mr')k(r', . (89) 

Working in the canonical ensemble in a finite volume V and inserting the 
expressions for p^"^^ and p, we get the remarkable identity 

/ H{r, r') dr' = {N - 1) p{r) + p{r) - Np{r) = . (90) 
Jv 

In other words, the integral operator H with kernel H{r,r') annihilates the con- 
stant functions and therefore cannot be inverted. The Ornstein-Zernike relation 
in its usual form ()87|) is thus unavailable in the canonical ensemble case. 



Acknowledgment 

We thank the two referees for a series of valuable remarks and suggestions. 



References 

[1] R. Evans: "Fluids in Model Pores or Cavities: The Influence of Confinement on 
Structure and Phase Behaviour", NATO Science Series C 529 (1999) 153 

[2] H. Loewen: "Density Functional Theory of Inhomogeneous Classical Fluids", 
J. Phys. Cond. Matter 14 (2002) 11897 

[3] J. D. Ramshaw: "Functional Derivative Relations for a Finite Non-Uniform Molec- 
ular Fluid in the Canonical Ensemble", Mol. Phys. 41 (1980) 219 

[4] J. A. White, A. Gonzalez, F. L. Roman, S. Velasco, "Density-Functional Theory 
of Inhomogeneous Fluids in the Canonical Ensemble", Phys. Rev. Lett. 84 (2000) 
1220 

[5] J. A. White, S. Velasco: "The Ornstein-Zernike Equation in the Canonical En- 
semble", Europhys. Lett. 54 (2001) 475 

[6] J. A. White, A. Gonzalez: "The Extended Variable Space Approach to Density 
Functional Theory in the Canonical Ensemble", J. Phys. Cond. Matter 14 (2002) 
11907 

[7] J. A. Hernando, L. Blum: "Density Functional Formalism in the Canonical En- 
semble", J. Phys. Cond. Matter 13 (2001) L577 



23 



[8] N. D. Mermin: "Thermal Properties of the Inhomogeneous Electron Gas", 
Phys. Rev. A 137 (1965) 1441 

[9] R. Lovett, C. Y. Mou, F. P. BufF: "The Structure of the Liquid- Vapor Interface", 
J. Chem. Phys. 65 (1976) 570 

[10] M. S. Wertheim: "Correlations in the Liquid- Vapor Interface", J. Chem. Phys. 65 
(1976) 2377 

[11] M. Requardt, H. J. Wagner: "(Infinite) Boundary Corrections for the LMBW- 
Equations and the TZ-Formula of Surface Tension in the Presence of Spontaneous 
Symmetry Breaking", Physica A 154 (1988) 183 

[12] M. Requardt, H. J. Wagner: "Does the Three-Dimensional Capillary Wave Model 

Lead to a Universally Valid and Pathology-Free Description of the Liquid- Vapor 
Interface near g = 07 A Controversial Point of View", J. Stat. Phys. 64 (1991) 807 

[13] M. Requardt, H. J. Wagner: "Implications of a Non- Analytic Small-fe Asymptoties 
of the Direct Correlation Function in the Liquid- Vapor Interface" , Physica A 196 
(1993) 59 

[14] R. Lovett: private communication 

[15] R. Lovett, F. P. BufF: "Examples of the Construction of Integral Equations in 
Equilibrium Statistical Mechanics from Invariance Principles", Physica A 172 
(1991) 147 

[16] J. L. Lebowitz: "Asymptotic Value of the Pair Distribution Near a Wall", 
Phys. Fluids 3 (1960) 64 

[17] M. Reed, B. Simon: Methods of Modern Mathematical Physics I, Academic Press, 
New York - San Francisco - London, 1972 

[18] P. Schofield, J. R. Henderson: "Statistical Mechanics of Inhomogeneous Fluids", 
Proc. Roy. Soc. London A 379 (1982) 231 

[19] J. R. Ray: "Correct Boltzmann Counting", Eur. J. Phys. 5 (1984) 219 
[20] T. L. Hill: Statistical Mechanics , Dover Publications, New York, 1987 



24 



